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I. INTRODUCTION 


Point-like interactions have provided a remarkably usefnl idealization for many situations 
in physics. In the context of scattering theory the concept of a point-like scatterer was 




introduced in 1934 by Bethe and Peierls l|. Fermi P| used and rehned their results to 
describe the motion of neutrons in hydrogenated substances (such as paraffin) by introducing 
what is now known as the ‘Fermi pseudopotential’. The idea is that when the scattering 
potential is concentrated on a very small scale rg (in the case studied by Fermi the range was 
that of nuclear interactions compared to the distances between the atoms) but nonetheless 
its influence on the motion cannot be neglected, one can characterize the scattering in a 
simple and efficient way by means of a few quantities like the scattering length, hnite in 
the ro ^ 0 limit. The problem of ‘how to separate the scales’ in the Schrodinger equation 
triggered by those 1930’s ^pers was addressed and elegantly solved over the years at different 
levels of formalism 3|, y, |5|, [y , [7| . The key to the solution relies in a proper dehnition of a 
‘delta-function interaction’ in dimensions greater than 1. 

Of course this paper will not be dealing with quantum mechanical scattering within 
matter, which is from many points of view a solved problem. The problem of ‘separation 
of scales’ however arises urgently in modern quantum field theory if stated as: “What is 
the quantum held theory response on length-scales L to a disturbance concentrated on a 
length scale rg -C T?” How does quantum held theory respond to topological defects and 
singularities, in particular of the metric? Once formalized in proper mathematical terms the 
two problems look much closer than one would think. 

The Casimir ehect falls in this class of problems. The penetration length (ro) of the 
electromagnetic held inside conductors is much smaller than the distance (L) between the 
conductors, which sets the scale of the experimentally measured force. We are interested 
in studying the dynamics of the conductors on the larger scale L by integrating out the 
electromagnetic held. This motivates the nomenclature ‘Casimir ehect’ for a much wider set 
of problems than Casimir’s original one. 

Another example: in any candidate theory of the quantum geometry of space-time the 
problem of dealing with point-like singularities will inevitably arise. Remember for example 
that the Ricci scalar for a point-like particle (like Schwartzschild’s solution) is a delta- 
function centered on the position of the particle. Quadratic huctuations of a non minimally- 
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coupled scalar (or of the metric) have hence a delta-function term in their Lagrangian. The 
effect of such a term must be considered together with the other known effects of the black 
hole metric. It is then of paramount importance to analyze the problem of how one or more 
concentrated singularities influence the spectrum and low-energy behavior of the fluctuations 
of the held. 

Analogous problems arise in condensed matter, quantum held theory and string theory 
since localized disturbances appear in all these theories, essentially only the names given to 
them being diherent (defects, domain walls, concentrated Aharonov-Bohm huxes and branes 
to name some). With this in mind we will set-up the problem in very general terms and 
even though not all the details map one-to-one on specihc examples the main results will 
apply to a wide class of examples. 

The main results of this paper are two. First, I will show that quantum huctuations of 
a scalar held 0 generate attractive forces between localized defects in the very same way 
Casimir forces act between metallic bodies. I calculate this force for an arbitrary number 
of defects with co-dimension^ 1, 2 and 3 (see Table I). Note that previously the Casimir 
ehect has been analyzed only for co-dimension 1. The main result of this paper is Eq. (I22D, 
which gives the interaction energy as a function of the scattering lengths of the defects and 
their relative separations. For co-dimension d > 4 the force disappears, as required by the 
properties of the self-adjoint extensions of the Laplace operator on the punctured (see 
3, Chap. X). 

Second, in the presence of two or more of these defects (of co-dimension > 1 ) the vacuum 
0 = 0 is hopelessly unstable and a localized tachyon mode is formed when the defects 
approach closer than a critical distance. At this critical distance the attractive force diverges. 
I calculate the wave function of the tachyon and show that it leads to condensation of the 
bulk held 0 to a vacuum expectation value 0 (x) 7 ^ 0 —but only in a limited region of space. 
The consequences of these observations for some models will be discussed in Section IVII 

^ Co-dimension is the number of dimensions transverse to a manifold. For example a point in 2 dimensions, 
a line in 3 dimensions and a surface in 4 dimensions all have co-dimension 2. 
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TABLE I: Flat manifolds divided according to dimension and co-dimension. The first line is the 
well-known Casimir problem, from the fourth line down the perturbation is ‘invisible’ to fluctua¬ 
tions. In this paper we will be dealing with manifolds in lines 2 and 3. 


II. THE INTERACTION ENERGY 

In this Section we will calculate the effective action of a scalar field coupled quadrati- 
cally to a static configuration of defects. The effective action Ses and Casimir energy S are 
proportional to each other 

SeS = -TS, ( 1 ) 

where T is the interaction time. In the following we will be interested in the Casimir energy 
of the problem. We will see that the part of the Casimir energy responsible for the interaction 
between the defects is a cutoff-independent quantity, meaning that the separation of scales 
can be performed effectively in this quantum held theory. 

We will consider the following action for the scalar held in d -|- 1 dimensions {h = c = 1): 



Here 6 is the d—dimensional Dirac’s delta-function, mimicking the concentrated disturbance 
on the held, /ij are constants, meaning that they do not depend on the held 0, but in general 
they can depend functionally on other helds living on the defect.^ The methods of will be 
used in order to dehne these d’s correctly. In this section, in order to keep things simple we 
restrict our attention to points in 1, 2 and 3 dimensions. We will add an arbitrary number 
of hat directions in Section El hence fulhlling our promise of studying co-dimensions 1, 2 


^ The /.ti’s are also called ‘brane tensions’. 
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and 3. 


Actions like Eq. © arise in different contexts. For example consider the case of a scalar 
bnlk field 0 conpled with N branes in cnrved or flat space (cnrvatnre can be easily inclnded 
in Eq. Q); or the case of a cosmic string (again the cnrvatnre ontside the string mnst 
be considered); or the case of electrons conpled with Aharonov-Bohm flnxes (in this case 
one has fermions rather than bosons bnt, after sqnaring the Dirac eqnation, the analysis 
is analogons [^). All these examples can be studied with the formalism introduced in this 
paper, so in full generality we will study the quantum fluctuations of the action Eq. ©• 
Some features of the Eq. (121) with only one delta function and constant m have been 

nnn 

studied before, for example in connection with cosmic strings scenarios The action 

© in one dimension with a single delta and a space-dependent mass term = m^{x) has 
been studied in In this paper we will consider the situation where a generic number 

of defects are present in d > 1 and is a constant. We will see that the situation will be 
different from that depicted in Ref.s 0, fl Q and unexpected physics is found. Moreover 
the rene ralization to x-dependent can be easily achieved by means of the techniques of 
Ref. [l^ and we will not comment on it in this paper. The inclusion of any other term in the 
action (|21) describing the dynamics of the surface itself would not affect our calculations. We 
consider the positions of the defects a* fixed and obtain an effective action. In the usual way 
this action can be used to describe adiabatically moving ai{t) (he. if the velocities \aij \ -C c). 

The forces between the defects can be calculated by taking the derivatives of the Casimir 
energy S with respect to its arguments {a*}. It will turn out that in general the forces are 
not additive, he. S is not a superposition of terms depending only on the relative distances 



We use the following integral representation of the zero-point energy of the scalar held 



( 3 ) 


where A is a cutoff and p is the spectral density of the Hamiltonian operator H 


N 



( 4 ) 


i=l 


^ The scale of variation of rn? should be of the order of the ‘long scale’ L. The purpose of this paper is to 
integrate out the physics at momenta ^ 1/ro which is symbolized by the delta functions in Eq. l(^. 
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Define also the unperturbed hamiltonian Ho as 


Ho = -V" + m\ 


( 5 ) 


Equation will look more familiar if the replacement E —> is done, and Ws are restored. 
One can obtain the spectral density p{E) as a functional of the propagator G{E) = 1/ {H—E) 
as 


p{E) = — lim ImTr G(E + ie). 

TV e^0+ 

In the following we will often write E + for E + ie when e —O’*". 

The propagator, Q{x',x;E) = {x'\G{E)\x) satishes the Schrodinger equation 

N 

(— Pi6{x' — Oj) — E)Q{x', x; E) = 6{x' — x), 

i=l 


( 6 ) 


( 7 ) 


and Qo = {x'\Go{E)\x) satishes the analogous equation without (5’s on the left-hand side. 
For m constant (which we will assume unless explicitly stated) and \m.E > 0 we have'^ 


Go{x\x]E) = < 


2VE-r 
rWi 


p 2 \/ E—m? 


if d = 1 


{H^o E - w?\x' - x\) ifd = 2 


( 8 ) 


E — rn^ \ x' —x\ 

Att\x'— x\ 


if d = 3, 


where H^'* is Hankel’s function of hrst kind of order 0. 


For d = 1 the problem is that of a scalar held on the line M in the background of a stack of 
d-functions centered on x = {oj} 1^. If we assume /i* > 0 they will attract each other, like 
metallic plates do via the Casimir ehect. These forces are not conhning and no new physics 
is obtained with the generalization obtained by adding n transverse directions. This is the 
usual Casimir problem. We will see how the situation changes dramatically when d > 1. 


To see how the solution for G is obtained, consider hrst the case with a single delta 
with strength pi = /r, placed at x = a. By solving the Lippman-Schwinger equation 
one hnds 


unction 

3 


la 


G{x , x; E) = Gq{,x , x; E) + 


1 


d — Go{cii Cl] E) 


Go{x\a]E)Go{a,x]E), 


( 9 ) 


^ In the following we will not use any special notation for vectors and we will indicate with \x\ the norm of 
a vector in 1,2 and 3 dimensions. 
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where a = —l//i. This solution is perfectly good in d = 1 and was the basis of the analysis 


in 


12j | for non-constant m?. For d > 1 it however suffers from a serious problem since 


^(a, a + r]I^ 
divergence j5[ 


—> cx) when the point splitting regulator r = |r| 
in a redehnition of a to obtain a hnite result 


0. One can reabsorb this 


Q{x\ X] E) = Qq{x\ x] E) -f 


ar-B(^\E) 


go(x\a; E)gii(a,x;E), 


( 10 ) 


where 




jy/E— 


if ^ 1 

if d = 3, 


( 11 ) 


47r 


where for d = 2 it has been necessary to introduce an arbitrary mass scale M which stays 
hnite when the point splitting regulator r —0. So a must be redehned such that when 
r —>• 0 

f 

a if d = 1 


— \ a -|- tT In Mr if d = 2 

ZTT 


( 12 ) 


a — 


47rr 


if d = 3, 


so the ‘renormalized’ a^. stays hnite. It is clear from this equation that one needs a positive 
divergent a to reabsorb the negative divergences when r —>■ 0. Large positive a mean 
negative very small fv (since a = —l//i). A small negative /i corresponds to a weakly 
attractive potential. Hence the point-like scatterer limit can be thought of as the limit of a 
concentrated attractive potential, zero outside a ^here of radius tq, with at most one bound 


state whose energy stays hnite when tq —> 0 


1 a ^JK 

m 


For d = 2 any attractive potential has 


at least a bound state and so we always hnd a bound state also for rn —> 0 (for d = 2 the 


flflQ); 


for d = 3 


dependence of a on M is reminiscent of a renormalization group how 
the bound state can be real or ‘virtual’ (he. a pole of the propagator G{E) located on the 
second Riemann sheet) its energy being hnite in the limit tq —> 0. The scattering length is 
(both for d = 2 and 3) a function of a and is hence hnite in the rg —0 limit. 

Another interpretation of this results comes from the theory of self-adjoint extensions of 


symmetric operators 


Here the ‘renormalized’ corresponds to a choice of self-adjoint 


extension for the Laplacian operator —A on the punctured 


In 


the self-adjoint 


extensions are not positive dehnite, meaning that they all have at least one (but it turns 
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out there is only one) negative eigenvalue. This corresponds to the bound state described in 
the paragraph above. In the punctured the self-adjoint extensions of —A can be either 
positive semi-dehnite (with a virtual state on the second Riemann sheet) or not (due to the 
existence of a single real and negative eigenvalue). 

The propagator with N deltas at positions {oj}, i = 1, and 1 < d < 3 can be found 


N 


g{x', X] E) = Qq{x, x]E) + '^ (T tti; E)go{aj,x] E). 

*,j=i 

The matrix T is dehned as (from now on we drop the subscript r on the a^’s) 


Tp- = (tti - B^'^\E)) 6ij - ^o(ai, aj] E), 


where 


Goiflii Ojj 1 E') 


0 


iii=j 


(13) 


(14) 


(16) 


It is now possible to explain why we limited our discussion to d < 3. The reason is that 


looking at thn 
self-adjoint 


aplacian A on the punctured one realizes that this operator is essentially 
, meaning that it has a unique self-adjoint extension: the trivial one. The 
4-dimensional delta function is ‘too small’ a perturbation to be seen by the Laplacian. What 
does go wrong in the renormalization procedure? The propagator in d = 4 is (|x' — x| = r) 


g^{x\x-E) = r) (16) 

oTrr 

~ — {^ E — -m^hir + O (1) ) if r ~ 0 (17) 

Stt / 

so we cannot choose a in an i?-independent way (because of the E — m? In r term) to 
remove completely the divergences as we did before. The low-energy limit of a 4—dimensional 
concentrated potential is hence trivial and we will not discuss this problem anymore. 

Having solved for the propagator we can hnd the density of states p simply by taking the 
trace and the imaginary part. The result® is 121 

(19) 


p{E) = pq{E) - -Im^lndetr(E-FiO+,{a}), 
TT OE 


® The only algebraic identity worth of notice is the fact that 

i:(r)5‘^6„(a.,a,;B) = -ATr l„r. 


( 18 ) 


'^3 
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where po = vr“^ImTr Go and the determinant of F is simply the determinant over the matrix 
indices ij. 

The term po in (jl9|l is independent of the presence, strengths at and relative positions of 
the delta fnnctions and we will neglect it in the following. The second term on the right-hand 
side of Eq. ( 111 , can be nsed to calcnlate the Casimir energy as a fnnction of the positions 
and strengths a of the scatterers: 


S = -Im 

271 


fA Q 

dE\fE-— In det T{E + zO’*', {a}). 
1 oE 


The interaction part of this energy is obtained by snbtracting from Eq. the 


( 20 ) 


same 


qnantity calcnlated with all the Lj,- = |aj — a 


oo. In this limit F becomes diagonal 


(considering that ImE > 0) and the energy (01 becomes a snm of self-energies of isolated 
objects. The self-energy of an isolated brane contains all the usnal nltra-violet divergences of 
the Casimir energy and mnst be treated with care Q| . The sharp (5-fnnction limit r —> 0 and 
the strong potential limit (sometimes called the Dirichlet limit) Vq —>■ oo are problematic for 
the Casimir energy already in d = 1 and the results depend on the order in which the r —> 0, 
Vq oo and A —oo limits are performed.® To avoid this problem we will consider the 
situation in which the potential V{x) is hnite, sufficiently smooth and localized over a hnite 
distance 0 < r L, where L is the separation between the branes. Physically, this means 
that we assumed the cutoff A of the held theory to be much larger than any other momentum 
(square) scale, like Vq or 1/r^. The divergences that then arise in the Casimir energy are 
entirely local: they will not affect the interaction energy. The procedure of subtracting the 
self-energies as described above then leaves a well-dehned, hnite interaction energy between 
the branes. After this subtraction is performed one can take the appropriate limits for the 
potential r —0 and Vq —>■ oo and, provided the interaction energy remains hnite as we prove 
below, the resulting interaction energy is unique. 

After performing the subtraction of the self-energies, the interaction energy can be written 


as 


S 




dEVE 



det T(E + z0+, {a}) 
det F(E -|- z0+, oo) 


( 21 ) 


® In this paper we are adopting a cutoff regularization of the quantum field theory. Different regularization 
schemes (like the widely used zeta-function, for example) give different results for the divergences. The 
finite part of the energy is, however, the same. 
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We keep using 8 to indicate the interaction energy, confident that this will not generate any 
confusion, since we will no longer be interested in the total energy. The integrand in Eq. m 
falls exponentially fast on the semicircle |E| —cx) of the complex E plane^ which allows us 
to integrate by parts, Wick-rotate to the negative E axis^ and send the cutoff A —> cxo. We 
can moreover remove the Im because all the quantities are real and positive on the negative 
real E axis (since the propagator Qq is real and positive for E real and below the spectrum) 
except for yJ—E + f0+ = i\/E. 

This leads us to a final, compact expression for 8 


8 = 


An 


dE j detr(— E, {a}) 
y/E detr(— E, cx)) 


( 22 ) 


This is the main result of this paper and together with the definition of T, Eq. (d , can 
be used to calculate the interaction energy of point-like scatterers due to fluctuations of the 
field 0. In the rest of this paper we present several examples of the applications of this 
formula. 


III. EXAMPLES 


As a first example and a check for our result, Eq. (I22D, let us calculate the well-known 
interaction energy between two delta functions at distance L, in 1 dimension (we assume 

Qi = 0(2 = cy 0 ); 


8 


1 dE ^ 
— / In 

47ryo VE 



^—2L\/ E-\-m? 


(1 — 2a\/E + m?Y 


(23) 


This formula reproduces the usual results for the Casimir energy of two penetrable plates 
in 1 dimension |l6j |. 

As another example in d = 1 consider the case of three repulsive delta functions («* = — 1 
for i = 1, 2, 3). The interaction energy can be calculated with the ease with which one can 
take a determinant of a 3 by 3 matrix. The result is plotted in Fig. ^as a function of the 
position X of one of the three deltas while the other two are held fixed at a: = 0 and a; = 5. 


^ In particular it goes to zero like e on the ray E = |i?|e*®, tt > 6 > 0, where L — min \ai — aj\ 

for any d. 

® During the Wick rotation we do not pick any pole contribution on the positive imaginary semi-plane of 
the first Riemann sheet because the total Hamiltonian Eq. 0 is self-adjoint. 
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FIG. 1: The interaction energy, in arbitrary units, for three delta functions on the line as a function 
of the position x of one of them and m = 0. a = —1 for all three deltas, one delta is held fixed at 
X = 0, another at x = 5. 

The interaction energy is not additive: the interaction energy of N semi-penetrable plates 
does not split into a sum of N{N — l)/2 terms due to pairwise interactions. Rather, by 
expanding the logarithm a reflection expansion is obtained in the spirit of Ref. 171. 

Before calculating the interaction energy for 2 or more deltas in d > 1 it is instructive to 
look at the case of a single delta function centered in x = 0, to introduce some properties of 
the bound state of a single delta. Consider the case d = 3. The propagator is (for Im > 0) 


^i\x'—x\%/ E—m''^ 

Q{x,x; E) = , , , -^-h 




(24) 


47r|x' —x| ' a _ iVE^- nE 167r2|x'||x| 

There is evidently a pole at E = Eq such that \/Eq — rn^ = —iAna. For a < 0 this is a real 
bound state at Eq = im? — and the wave function 'ipQ of this bound state is obtained 

by noticing that 


E.-E 


xpl{x)'ipo{x), 


for E near the pole Eq. We hence expand (121 about Eq to hnd 

■M- 


'tpo{x) = 


-a) 


47r(—a)|x| 


\X\ 


(25) 


(26) 


For CK > 0, on the contrary, the pole is on the 2”'^ Riemann sheet and hence is a virtual state 
and does not belong to the spectrum of H. Whether this pole is real or virtual, physically 
it represents the s-wave scattering over a concentrated attractive potential, l/a is indeed 
proportional to the scattering length in the s-wave channel j^. The s-wave is the only 
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0.5 1 1.5 2 2.5 3 

L/Le 


FIG. 2: Interaction energy £ (the continuous line is Rei5 and the dashed line is —Imi5) in units of 
1/Lc, for two delta functions as a function of their distance L. a) The case with m = 0. b) The 
case with m/M = 2. 


contribution surviving in the limit when the scatterer is small compared to the wavelength 
l/VE-mh 

We have to require the spectrum of H to be contained in the positive real axis for the 
vacuum, 0 = 0, of our held theory to be stable. So if a < 0 we have to choose m > 47 r(—a). 
If a > 0 any choice of m, in particular m = 0, is enough to ensure the stability of the 0 = 0 
vacuum.® 

Considering the case d = 3 further, let us now calculate the interaction energy between 
two identical delta functions with ai = a 2 = a > 0 (so, according to the previous paragraph 
no bound state exists for isolated scatterers) at a distance |ai — 02 ] = L. After the Wick 
rotation and dehning k = y/E we obtain 
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(dvra + y/k'^ + m?) 


(27) 


It is not difficult to see that it exists a critical distance L^, being the positive solution of the 
equation 

r ^ _ I _ 

dvra + m’ 

such that if L < Lc the argument of the logarithm in Eq. (EZD becomes negative for suf- 
hciently small k and we get a negative imaginary part in the Casimir energy. A negative 



® It is worth noting that the opposite choice for the sign of a is needed to avoid a bound state in the d = 1 


case. 
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imaginary part of 8 means, as usual, an instability of the 0 = 0 vacuum in the presence 
of the two (5’s. In fact, by studying the eigenvalues of the matrix F one can see that for 
L < Lc the spectrum of H has a bound state with negative E, and since E = ui'^ this is a 
clear indication for the existence of a tachyon. We will return on the implications of this 
instability for the low-energy physics. 

The force E = —OE/dL, always attractive and central, diverges logarithmically at the 
critical length Lc- For m = 0 and {L — Lc)/Lc. <C 1 one hnds: 






In 


L-E 


(29) 


The long-distance behavior, L S> of the force depends on the mass of the boson 0. 
For m > 0 the potential between the two (5’s decreases exponentially. For m = 0, instead, a 
power-law tail is obtained: 

This 1/E law is stronger than the Casimir-Polder law (induced polarization interaction 
which falls like 1/L^. This means that we should not think of these delta functions 
as mimicking polarizable molecules or metallic particles. Indeed, to correctly describe a 
metallic sphere of radius i?, surface S and penetration depth tq one should rather assume 
that ro -C i? adding hence to Hq in Eq. a potential lF(x) = — = jji5{r — R) 

and send y, ^ oo before sending i? —0. This clearly is a different limit than the one we are 
describing here. 

Now that we have discussed the divergences associated with F, its renormalization and we 
know about the existence of vacuum instabilities related to negative E bound states of the 
hamiltonian El we are ready to tackle the two dimensional case where all these complications 
arise at the same time. 

The free propagator is Qq{x/x] E) = jHq^\EE — rrElx' — a;|) if \m.E > 0. Notice that 
even for a single delta no choice of a eliminates the bound state. There will always be at 
least one bound state with energy Eq = m? — This is due to the fact that any 

attractive potential in 2 dimensions has a bound state. We must choose our mass such that 
Eq > 0 and the instability is not present (it suffices that m > However we will see 

that in d = 2, exactly as in the d = 3 case discussed above, in the presence of two or more 
d’s there exists a critical distance such that for closer approach a bound state has E < 0, 
generating a tachyon again. 
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Take N point-like scatterers, each with renormalized strength and a renormalization 
mass M. It is convenient to dehne Mi = so T in m is 

~ IM- ) ~ Ojj] E). (31) 

The interaction energy for two identical deltas (Mi = M 2 = Ai, and m > Ai as required 
for the stability of isolated scatterers) separated by a distance L is {k = a/E) 

s^^rdkiJi- 

271 Jo \ In^ {Vk^ + mV-M) 

where Kq is a Bessel function K of order 0, and the critical length is the solution of the 
equation Ko{mLc) = Inm/Ai. The force diverges as L —in d = 2 as well but the 
explicit expression is more difficult to recover. For L ^ the force is exponentially small, 
since we had to assume a mass m > 0 for the held 0. 



IV. LOCALIZED VACUUM INSTABILITY 


Let us calculate the shape of the tachyon in 3 dimension found in the discussion after Eq. 
(fTfjl (take m = 0). Let us hrst notice that (for any number of scatterers) the Wick-rotated 
matrix r(— E) is real and symmetric and can hence be put in diagonal form. In the case 
at hand we have only two delta functions with equal strength a and one can show that the 
spectral decomposition of T”^ is 


where 


and 


(T i)ij(-E, {01,02}) = 

I1I2 


71 

72 


q-lVe 

AttL 

q-lVe 

AttL 


+ 


Ve 

dvr 
4 TT 


a. 


a, 


^(1) 

vi^) 


{—-—} 

'V2’ 


(33) 

(34) 
(36) 

(36) 

(37) 
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The bound state pole is generated by a zero E* in the 71 eigenvalue. For L < Lc = l/dTro; 
we have E* as the real positive solution of the equation 71 (i?) = 0 (remember: the integration 
variable E appears as —E in F so positive E here are real, negative eigenvalues of H) 


Ve- 


^—L'/E 


Lc L 

Comparing the behavior of the propagator for E close to a pole E* 


Q{x,x] E) ~ 


E* - E 


(38) 


(39) 


with 

g(x',x-,E) ~ — (40) 

we hnd the wave function of the (not normalized) bound state as 

g-FF*|x-ai| g-\dE*|x-a2| 

—i-i-1-i-i— 

\x — ai\ \x — a2\ 

Something can be said also in the case in which we have many identical defects (and 
assume all the Oj/s are of the same order of magnitude), without necessarily having to 
solve the equations explicitly. Instructed by the previous analysis we can state that the 
gronnd state will be a highly symmetric state v ~ { 1 /V^,..., 1 /V^} which will then give 
a symmetric wave fnnction '^o(a^) oc '^^ViQ{x,ai] E*) delocalized over the entire array of 
defects. The positivity of the nds coincide with the constraint that the gronnd state must 
not have any node. 

Hence for L < Lc & free field theory conpled to these defects does not make any sense. 
Its vacnnm state 0 = 0 is unstable. The imaginary part of the energy (as an analytic 
continnation to L < Lc) is related to the ‘decay time’ of the vacnnm state, due to particle 
creation. 

Let us for a moment speculate on the conseqnences of this instability. Adding higher 
order terms in 0 to the Lagrangian —one can for example think of adding a term— 

should eventually stabilize the held with a vacuum expectation value (vev) (j){x) 7 ^ 0 in a 
somewhat large region around the two scatterers. However the actnal valne of the vev (j){x) 
and the size and shape of the condensation region cannot be easily constrncted and will be 
subject of future work. 

This scenario of a local condensation and creation of localized vacnnm instabilities dne 
to defects could be interesting in inhation cosmology as well (the held 0 being the inhaton). 
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It must be also remarked that a similar scenario occurs in brane cosmology when an open 
string has its ends attached to two D-branes When the branes are pushed closer than 
a critical length one of the modes of the string becomes a tachyon. 


V. EXTENSION TO n - 1 TRANSVERSE DIMENSIONS 


Now that we know the density of states p{E) for the ‘basic’ problem of points in 1, 2 and 
3 dimensions, we can move along the lines of Table I to generate solutions for manifolds with 
co-dimensions 1, 2 and 3. We shall then add n — 1 transverse, flat dimensions. The total 
dimension of the space is now d + n — 1. The calculations in the preceding part of this paper 


201 and 


can be recovered by putting n = 1 in all the formulas. We will use the methods of 

where one solves for the density of states p{E) of the basic problem on a d—dimensional 
section and insert the result in the equation for the energy per unit n — 1-dimensional ‘area’ 

^ 0 

/ dE^(x/p^ + E-x/p^]p{E). (42) 


(27r)-i 7o 2 


The subtraction — removes a divergent but a—independent term, since the integral 
j dEp{E) is a-independent. We will also remove the a-independent ‘self-energy’ terms by 
subtracting from p{E) the density p{E, oo) with all aij —>• oo. We can then perform the 
(dimensionally regularized) integral over p, Wick-rotate and perform an integration by parts 
on E to obtain 


g{n) ^ 


1 r(i- 


f)sin^ 


detr(-E;{a}) 


/o det r(— i?; cx)) 


(43) 


27r (47r)’^/2 

For example the interaction energy (per unit length) of two straight, inhnite strings in 3 


dimensions (d = 2, n = 2 so d -1- n — 1 = 3) put at a distance L is 




Stt 


dE\n 1 - 


Kl{L^E + m?) 

-|- m ^/M) 


(44) 


One can hence calculate the interaction energy of any two flat manifolds due to the quantum 
fluctuations of a bulk held 0. As an example in co-dimension 1 consider the Randall-Sundrum 


scenario 


21| with two branes at a distance from each other. The huctuations of a given 


component of the metric G or of a bulk held </> (see 2^ and references therein), have a space- 
dependent mass with two delta functions singularities on the two branes. The attractive 
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force due to the quantum fluctuations of this held has a Casimir-like behavior. The curvature 
in the 5*^ direction does not change the physics. If however the branes have co-dimension 
2 or 3 (and are dehned as the limit of an attractive potential) is in the class of problems 
that we have studied in this paper and a perturbation would eventually condense, if Tc < L^. 
Following the same arguments above we can also say that if the branes have co-dimension 
> 4 the huctuations in the bulk will not see the brane. The cosmological implications of 
such a scenario will be subject of future work. 


VI. OMISSIONS AND APPLICATIONS 

The propagator, Eq. m, comes directly from scattering theory. In that context it was 
natural to assume that the interaction between the particle and the scatterer (consider 
Fermi P| and Zel’dovich P| examples) is attractive. One considers an attractive center 
whose attraction grows when ro —0 such that at most one bound state is present and its 
energy remains hnite {i.e. of (9 (1)). Even thought we assumed that only one bound state 
is present at energies of O {1) this is the most generic situation that can occur in scattering 
theory. In fact if a second bound state is present, it will be an energy O (1/rg) below our 
bound state. In the limit rg —0 its influence on low-energy scattering disappears. It goes 
out of the spectrum. In scattering theory however such a negative energy state is harmless. 
This is not the case for a bosonic held, for which it represents a tachyon. 

One may wonder: what happens if the potential is repulsive and concentrated? The 
answer is that for d > 1 its inhuence on the scattering matrix (and hence on the spectrum) 
disappears when tq 0. Obviously this is not true in 1 dimension because we cannot ‘go 
around’ the scatterer. For a repulsive potential in d > 1 the renormalization procedure 
leading to (USD cannot be performed since fv has the wrong sign and sending rg —0 just 
kills the correction to Qo in m- 

More precisely, if in the Lagrangian we include a term Vod(rg — \x\)(j)‘^{x) with Vq > 0 
and then we take the limit Vq —>■ cx) and rg —0 with Vorg hnite, the spectrum we obtain is 
just the free one: the scatterer disappears. In a sense, the only smile the Cheshire cat can 
leave behind is the lightly (he. O (1) instead of O (1/rg)) bound (or virtual) state. If this is 

If the brane is inside an horizon for the 5d metric most probably this assertion is not true, however. But 
this is not the case for the Randall-Sundrum model. 
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not present then the scatterer is invisible to the fluctuations.^^ If the purpose of calculating 
the effective action was to calculate quantum corrections to a classical solution (as often 
occurs) then we deduce that for a repulsive potential or for d > 4 the classical solutions are 
unchanged by quantum fluctuations. 

Let us now comment on two possible applications of the formalism we have developed: 
cosmic strings and concentrated Aharonov-Bohm fluxes. We anticipate that further work is 
required in both cases. In the literature on cosmic strings the difficulty generate d by a bound 


0 , 0 . 


In that 


state tied to a single cosmic string has been recognized a long time ago 
context the bound state arising from Eq. m is rightly considered fictitious, because the 
smoothed potential is always positive (/i > 0). Nonetheless, in [0 after projecting out this 
bound state at Eq < 0, the propagator m is trusted and shown to be in good agreement 
with the numerical solution of the smoothed problem. It is not clear if projecting out a 
state from the propagator by hand has non-trivial (wrong) consequences on the density of 
states and the Casimir energy so we preferred not to follow this path even if it gives correct 
results for other quantities. We hence required the field to have a non-zero mass so that this 
bound state is stable. In the end it is not clear if the Casimir attraction and the birth of 
the tachyon could arise in cosmic strings coupled with bulk fields. 

Another example to which the above techniques and results should be relevant is the case 
of a fermion around a concentrated tube of flux (Aharonov-Bohm case^^). The spectrum of 
Dirac’s equation can be inferred from that of a Klein-Gordon equation after squaring the 
former. The fact that we are dealing with fermions rather than bosons, is not a difficulty. 
Another difficulty however arises: for Aharonov-Bohm fluxes and the more^neral case 
of cosmic strings charged under some U(l) symmetry, it has been shown 2^ that the 
contribution to the scattering cross section given by the non-zero external vector potential 
is asymptotically larger in the low energy regime than the contribution of the singularity 
in the core. Since we believe that cross sections and Casimir forces are tightly bounded 
quantities we would not apply any of the above arguments without treating the propagation 


We have already remarked on the impossibility to consider the Casimir-Polder interaction between small 
metallic spheres as the interaction of point-like scatterers the way we construct them here. Small metallic 
particles of radius r still have a penetration depth ro r so effectively are co-dimension 1 surfaces. 

The one loop energy of QED flux tubes has been calculated in |2^ using a combination of analytical and 
numerical methods. Our method could be used to calculate the interaction energy of two such tubes in 
the limit where their radius is small compared to their relative distance. 


18 




in the external space properly. This will be done elsewhere. 

Renormalization of branes coupling for a single brane (or ^-function in our case) with 
co-dimension 2 (and dimension 5) in a conical space has been studied in Arising from 

local divergences, the renormalization flow is not affected by the presence of other branes and 


in Q 


the results in [la apply also to our situation. Their renormalization of the brane coupling 
H (A 2 in their notation) corresponds to our renormalization a —> ar- Their renormalization 
of the effective action corresponds to our subtraction of the a-independent terms in the 
Casimir energy. These two are the only subtractions needed (if 0^ terms are not present) 


and it is heartening to see that our results coincide with those of [l^. Moreover one can 
make an amusing observation if one compares the two approaches to the delta function, the 
one in terms of scattering (that we used here) and the one in 1^ in terms of renormalization 
group. Notice that the renormalization group flow for p is IR free and has a Landau pole: the 
location of the Landau pole coincides with the location of the bound state in our approach. 


VII. CONCLUSIONS 

We have calculated the force between an arbitrary number of surfaces (branes) with co¬ 
dimension > 1 due to the quadratic fluctuations of a boson 0 living in the bulk. The force 
turns out to be attractive and it diverges when the distance between the branes approaches 
a critical value Lc- This phenomenon has no analogues in the widely studied co-dimension 
1 case. 

The divergence of the force is accompanied by the birth of a vacuum instability, a mode 
with negative mass squared localized around the scatterer. In 3 dimensions, the long-range 
properties of this force (decreasing like 1/L®) are shown to be different from the Casimir- 
Polder 1/L^ law, the explanation relying in the proper mathematical dehnition of the point- 
like limit. 

Some implications of these effects have been pointed out. 

See also 0. 
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